arXiv:hep-ph/9701378vl 24 Jan 1997 


UM-TH-97-05 


|hep-ph/970rFf8 

February 1, 2008 


Renormalon Variety in Deep Inelastic Scattering. 


R. Akhoury and V.I. Zakharov 

The Randall Laboratory of Physics 
University of Michigan 
Ann Arbor, MI 48109-1120 

Abstract 

We discuss the renormalon-based approach to power corrections in 
non-singlet deep inelastic scattering structure functions and compare it with 
the general operator product expansion. The renormalon technique and its 
variations relate the power corrections directly to infrared-sensitive parame¬ 
ters such as the position of the Landau pole Aqcd or the infinitesimal gluon 
mass A. In terms of the standard OPE these techniques unify evaluations 
of the coefficient functions and of matrix elements. We argue that in case 
of deep inelastic scattering there is a proliferation of competeing infrared 
sensitive parameters. In particular we consider the gluon and quark masses, 
virtuality of quarks and Aqcd as possible infrared cut offs and compare the 
emerging results. In the standard renormalon technique where Aqcd is the 
infrared parameter, the argument of the running coupling is crucial to ob¬ 
tain the correct x dependance of the structure functions. Finally we discuss 
the limitations of the use of the renormalon based methods for determining 
of the X dependance of the power corrections. 
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1 Introduction 


Generically, renormalons allow for a parametrization of infrared sensitive contributions to 
hard processes in QCD within the framework of perturbation theory (for the basic ideas see 
iBi and references therein ). If one considers only infrared safe observables then the 
infrared sensitive contributions are power suppressed. Moreover, in many cases it does not 
matter which particular infrared parameter is used in calculations. It can be, for example, 
the position of the Landau pole in the running coupling, Aqcd, or (at the one-loop level) 
a hctitious gluon mass A. Since renormalons are a pure perturbative construct the answer 
is obtained directly in terms of the infrared parameters chosen. If, on the other hand, the 
treatment of the same process is possible within an operator product expansion the evaluation 
of the power corrections is a two-stage procedure. First, one calculates perturbatively a 
coefficient function in front of operators. The matrix elements of these operators provide 
then a measure of the infrared contribution, both perturbative and non-perturbative. Thus, 
renormalons £x the matrix elements perturbatively. Finally, one adjusts the overall scale of 
the matrix element to allow for its non-perturbative enhancement. This tacit assumption is 
behind almost all the applications of renormalons Q. 

The best known example of this kind is the gluon condensate (0 |q!s(G“^)^| 0). It can be 


treated non-perturbatively [^. On the other hand if one introduces a hnite gluon mass A 
then 0 : 



where we keep only the term nonanalytic in A^ since only such terms can be consistently 


attributed to the infrared region. Instead of introducing A 7 ^ 0 one can evaluate the gluon 
condensate associated with a renormalon chain 



ren 



( 2 ) 


where asik"^) is the running coupling corresponding to a single-term /3-function and the 
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integration over the pole in as{k‘^) has been defined as the principal value (p.v.). In this 
example the use of renormalons, i.e. of Eqs. (i, 0. does not enhance the predictive power 
of the theory. Indeed, Eqs. (0,(0) cannot be taken literally and one reserves for an unknown 
rescaling of the matrix element to fit the data. What renormalons do achieve is a short cut 
of the procedure since, say, terms of order XHnX^ can be evaluated directly for physical 
quantities without invoking the operator product expansion as an intermediate step. This 
may be a decisive advantage of renormalons in cases where there is no OPE available. 

Thus at first sight the alternative is that either there is an OPE and then the renormalons 
are a particular model for the matrix elements involved or there is no OPE and then the 
renormalons could be a substitution for it. This logic appears to be defied by the example 
of deep inelastic scattering. Indeed, the power corrections can be treated either via the 
OPE [0, 0,1^ and or via renormalons [1^, |T0]. While the OPE reserves for at least 

one unknown matrix element for each moment, renormalons are claimed to fix the whole x 
dependence of the 1/Q^ correction Moreover the same approach can be used to 

study the power corrections to fragmentation functions . 

Motivated by these very interesting observations we look, in this paper JT^, into the 
anatomy of the infrared sensitive contributions in the renormalon calculus as applied to DIS. 
We find that already at the one-loop level there are a variety of possible infrared procedures 
which do not obviously give the same answer. For example, as discussed in more detail later 
in the paper, one possible strategy underlying the applications of renormalons is to reduce the 
matrix elements governing higher twist contributions to those of the leading twist. Typically, 
the reduction produces a factor / proportional to: 

rl 


/(e^m^A^) ~ [ dyX{y)ln{X{y)) 

JO 


where 


771 ^ X^ 

x{y) = + + 

A, m are the gluon and quark masses respectively and e is the virtuality of the quark, = 
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e^. This proliferation of the infrared sensitive parameters, that is A,m,e, is due to the 
fact that there are both soft gluon and soft quark lines. On the other hand the classical 
applications of renormalons assume that there is only a single gluon line made soft through 
an insertion of a renormalon chain, while the other lines are hard. Thus in case of DIS, the 
renormalon-related parameter A determines the result only upon forcing the other parameters 
(m, e) to vanish - m^, -C A^, which is not necessarily natural. One could of course assume, 
for example, that to the contrary, S> A^. In addition to A and m, one could also 
explore the infrared sensitivity by the standard renormalon chain technique in which the 
power corrections are proportional to ^qcd^/Q"^■, etc. In sections 2 , 3 and 4 we consider i 


m 


more detail the x dependance of the non-singlet structure functions with A, Aqcd and m , 
respectively, as dominating infrared parameters. Section 2 essentially reproduces the results 
of 0,0 0. In section 3 we hnd that not only is the predicted x dependance a function 
of what the argument of the running is, but in addition, two arbitrary scales typihed by two 
different integrals multiplying Aqcd^/Q"^ are introduced, in general. This difference in fact 
arises because some contributions come from Feynman integrals which are collinear divergent 
whereas others from ones that are not. This is to be contrasted with the case discussed in 
section 2 where there is just one unknown parameter A and is an indication of an inherent 
infrared instability. As for the dependence on the virtuality e it also indicates that in fact 
we are dealing with an infrared unsafe quantity. This would, in general, become manifest 
in the two-loop approximation for the power corrections. To see this one does not need to 
evaluate the second loop explicitly but it is enough to note that the anomalous dimensions 
of the operators governing the leading and higher twist contributions are different. This is 
discussed in section 4 . In section 5 we present our conclusions. 
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2 Gluon Mass as an Infrared Cut- Off. 


As is mentioned in the introduction there exist a variety of choices of infrared sensitive pa¬ 


rameters. Thus far the dispersive approach to the running coupling [|^, and the renormalon 
chain in the large Nf limit have been tried , with identical conclusions for the x depen- 
dance of the DIS structure functions wherever there is overlap. To these we will add the 
cases of A, m 7 ^ 0 and of the Landau pole contribution in the running coupling. These cases 
are considered in this and subsequent sections. 

To introduce notations, we consider the partonic structure tensor dehned as 


nUp.i) = (3) 

where is the amplitude for + q —>■ q' + g + ... where q is the parton ( quark for our 
purposes ) of momentum p. An average over the initial quark spins is understood and the 
standard decomposition of the is used: 








Fl{z,Q^) + 


( 4 ) 




p-q 


+ Puq^,) + Pm. 


{p-qf\ F,{z,Q^) 


^ 2 
qz 


p-q 


where z = j‘Ip ■ q, q^ = —Q^ and = 0. Note that we put the mass of the quark 
exactly zero while keeping a hnite gluon mass so as to ensure the role of A 7 ^ 0 as a unique 
infrared cut off. For the structure function Fl only the diagram with emission of a real gluon 
contributes to lowest order in For single gluon emission we have : 


J -^^{2Try6+{{k + g)^)5+((p - kf - 

( 5 ) 

where in the Feynman gauge dp^- = —Ppa- 

To perform the integral we use a Sudakov parametrization: 

I p2 

k^ = np^ -I--— -n^ -I- k^ (6) 

2k 
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where = 0, = 0, n • /cj_ = 0, n ■ p = 1, p ■ fc_L = 0 and we may choose 


= 


+ zp'^ 
p-q 


( 7 ) 


Moreover, for these variables: 


(Tk = 


TT r dn 


dk^dki 


K 


( 8 ) 


For detecting the non-analytical terms the limits of the k"^ integration are also important: 

X^z 

- < fc" < 

2 ; 


1-z 

It is straightforward to hnd for F^: 


< e < 91- 


(9) 


Fr = 2C, 


as 4^^ 
dvr 


n2 \2\2 A k"^ ~ 




( 10 ) 


with the limits of integration as indicated in (^. From (|10D we immediately get for the term 
independent of A and representing therefore the leading twist: 

(y.Q 


Fl = Cf-^Az 
dvr 

which is a well known result. Furthermore for the terms X^lnX^ and XHnX^ we get: 

a. Az'^ 


(F. 


L)\2 


= c 


and 


{FL)x‘i — Cp 


'2tt g2 
a. Az^ , 


■ 2X^lnX^ 


-?>)XHnX^ 


( 11 ) 


( 12 ) 


(13) 


T3|) are understood in the 


27r 

respectively. It is worth emphasizing again that equations 
sense that only non-analytic terms in are kept. 

The results (1^ and (p!3D do reproduce the predictions based on the dispersive approach 


to the running coupling []n|. This coincidence of the results comes as no surprise since the 
dispersive formulation is, in fact, aimed at improving the high-energy behaviour of the theory 
with massive gluons. However, the change in this high-energy behaviour obviously does not 
affect the infrared sensitive pieces which are non-analytic in A^. 
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3 Landau-Pole Parametrization. 


As the next way to parametrize the infrared sensitive contributions we examine in this 
section the contribution of the Landau pole in the running coupling, in which case the power 
corrections are proportional to {Kqc■ For dehniteness we again consider hrst the 
longitudinal structure function. Since now the infrared sensitive parameter will be Aqcd 
we can put A = 0. Moreover, we now account for the running of the effective coupling 
while still using the kinematics of one-gluon emission. The crucial issue then is what is the 
argument of the running coupling. We perform hrst the integration assuming that it is 
that determines the running, i.e, as{k‘j_). We shall see however that to reproduce the results 
of the previous sections one should assume that the true argument of the effective coupling 
is kl/{l-z). 

To depict the contribution to of the Landau pole we hrst integrate over k and k"^ (see 


Eqs. 


Thus we start with 

42;2 /. 




0^s{k±) ^^2 f ~ 


dvr 


K k^ 


K 


6ie + -Q^ 






K 


K 


K / 

' e + kj ^ 

K Z 


(14) 


Then the integration over k"^ is trivial while the integration over k can be expressed in terms 
of K±. 

2n± = 1 + z ± Y^(l — zy — 4{k‘^/Q‘^)z{l — z). (15) 

In this way we come to the following integral over k‘j_\ 

4^2 rQHl-z)/4z a,{kl) dkl 


Fr = 2 C^ 


F ■ 


Q'^ Jo 


dvr 


(16) 


'il-zy-4:{kl/Q^)z{l-z) 

In order to pick up the renormalon contributions we substitute the following representa¬ 
tion in Eq. (|I^ : 

/ 1-2 \-^dl 

(17) 


asik]) = da 


J^QCD, 
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where j3i is the hrst coefficient in the (3 — function. Furthermore, to perform the integration 
over k± it is convenient to introduce a new variable, y = 4^(1 — which allows us 

to disentangle the z and kx_ dependences: 

^2 / A 2 \ O'/^l 


Fr, = 


2Cf 


An 1 — z 

Finally we get for Fl the following expression 
2Cf Az^ 




fl - zV-F r(l - aA)F(l/2) 


Fr = 


An 


z Jo 




Az 


F(3/2 - a/3i) 


(19) 


which exibits the infrared renormalons corresponding to the poles of F(1 — Since 

~ exp{—a/as{Q‘^/AQQjj)) the corresponding power ambiguities are propor¬ 
tional to (Aqcd/<5^)> {AgcD/Qy- 

In order to extract the power behaviour one must dehne the integral over a using some 
prescription. We next outline how this integral may be dehned via, say, a principal value 
prescription for the poles at a = n/Pi, n = 1,2,.... First we note that small values of a 
correspond to large values of k± and hence for this we get the usual renormalization group 
improved perturbative answer . Keeping this in mind we may divide the integration region 
thus: 

Q r? — Q T) -\- Q 

( 20 ) 


a e |0 ,«d| = |0,^] + i:i-;-1 

W n>l W 


n>l /^1 

where, 0 < s < 1, and a typical choice could be s = 1/2. The hrst region gives the 
perturbative answer, as we have checked, and the rest give the power corrections in the 
principal value prescription. To make this explicit, consider the contribution from the pole 


at a = n/f3i to the integral (see Eq. (|T^): 


da 


r(l-aA)r(l/2) 


Xl-z)Q\ 


F(3/2 - cr/3i) 


which is dehned by: 


n + s 


^zAlcp \ F(l-aA)F(l/2) 

F(3/2 - uA) 


isf 


( 21 ) 


( 22 ) 


















This can be written after some algebra as: 


Xl-z)Q^) - 


(23) 


where, the integral is defined by : 

ir = (- 1 )"^ 


From this, for the first two leading power corrections we hnally get: 

Cf 4 ^^ ^%CD 


/■* da 



Jo simra 

r(l + n + cT)r(3/2 — n — a) 

r(l + n — (T)r(3/2 — n + a) 


(24) 


(Al)i/Q 2 - 


(Al)i/Q4 - 


27r 1 — z 
Cf 16^3 


’ jpv 


27r (1 


QCD jpv 

" ■ 


(26) 


We would like to emphasize that we could have chosen some other prescription for dehning 
the integral over a and hence the overall scale of the power corrections is arbitrary. An 


important aspect of Eq. (^) is that the z-dependence of the power corrections differs from 
the corresponding z-dependence of the power corrections due to A 7 ^ 0 at the one-loop level 
(see Eqs. (|I^ , (|T^). Thus, the predictions for the ^-dependence of the power corrections 
are sensitive to the exact argument of the running coupling . 

Let us focus first on the question concerning the most sensitive point of the derivation. 
Imagine that the true argument of the running coupling is in fact k\/{l — z), i.e,: 


as{k^_) 


a 




(26) 


This would correspond to a change in Aqcd in the Eq. (0): 


^^QCD 


{l-z)-A^ 


QCD 


(27) 


which clearly would bring Eq.(p5|) in line with the predictions based on 7 ^ 0, see Eqs. 
(0, (^). The corresponding principal value integral will be denoted by and for future 
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reference it is given by: 


ir = (-1)^ 


. TT' 


3/2 


/■* da 

(4z+ <t) 

a) 1 

Jo sirnia 

r(l + n + cr)r(3/2 — n — a) 

r(l + n — (T)r(3/2 — n + cr) 


(28) 


Since Fi is not logarithmically enhanced in the leading order, one may think that there is 
no reason, a priori, to believe that the coupling in Eq. ([T6|) should run like aa{k\). However, 
since it is well known that for the structure function F 2 , perturbation theory upto the 
two loop order for the leading contributions which are logarithmically enhanced as z —>■ 1 is 
consistent with the identihcation let us consider the 1/Q^ corrections to it. Here we 

will hnd not only the pattern indicated above but in addition a different kind of arbritariness 
related to the fact that the structure function F 2 recieves contributions which are collinear 
divergent. 

Using dimensional regularization to isolate the collinear divergences we readily obtain 
the following expression for the real contribution to F 2 : 

, _ 2 e 1 


preal ^ 2C 


-k 


27r h ^{i_zY_4^{kl/Q2)z{l-z) 


Ota 




,Qyi-£) ,, 

4z dk I 


-2e 


-C, 


as 3z 1 

2711 — z Q'^ Jo 


42 


dk\ 




4z 


dk I 


;_1_ 

\/(l - - 4:{kl/Q^)z{l - z) 

1 

il-zy-Aikl/Q^)z{l-z)' 


(29) 


There are of course the virtual contributions which in particular tranform factors like 
into ”+ ” distributions, however, we omit them for the moment as they are not crucial to 
the argument. It is easy to check that integration over k± reproduces the well known result 
for this structure function to leading twist. Note that only the integral multipying the hrst 
two terms produces a collinear divergence proportional to 1/e ( e is related to the number of 
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space time dimensions {D) hjD = 4 — 2e) whereas the second type of integral does not. 

In order to determine the Aqcd'^/Q'^ contributions from the above using the renormalon 
method, let us hrst consider the result when we take as to be as(Y^) for all the terms on 
the right hand side of Eq.(|29|). Proceeding as for the case of Fl we hnd: 


preal ^F 

^ “ 27r 




(30) 


where J and I are two different types of integrals given by: 


^ r(—(T/3 i — e)r(l/2) 
r(l /2 - aA - e) 

r(i - crA)r(i/ 2 ) 


J = 


I = 


da 


da 


J 




(31) 

(32) 


J r(3/2 - a/3i) • 

Notice the collinear divergence in J for a ~ 0 which was identihed with the perturbative 
region. One may imagine dehning J and I by two different prescriptions for obtaining the 
power corrections ( for which the ratios J/I would be different ) and in this way we see 
that two arbitrary scales would appear in the predictions. Thus if and denote the 
appropriately dehned integrals in the integration region over a which produces the power 
correction (1/Q^)'^, then: 


Tpreal 

^2 


1/Q2 


Cf 

2ti 


1-z 


— (1 + ^) j 4^Ti + ( — 




A, 


QCD 


In the principal value prescription, for example, is different from (see Eq. 
is given by : 


(33) 

and 


Jr = (-1)' 


. TT' 


3/2 


da 


Pi Jo simra 


(4z- 


^QCD 

Q2 


y{n-l/2 + a) 


(4z- 


■QCD \-a, 

Q2 ) ' 


n 


1/2 


a 


r(l + n + a)T{?>/2 — n — a) r(l + n — a)T{?>/2 — n + a) 

(3 

If we would have used the running as{k‘]P instead of as{j^) then in particular, we would 
have an extra factor of 1 — z in the denominators of Eq. (|33|). The z dependance of the 


higher twist contributions as given in Eq. (|3^) is of a similar type as that obtained from 
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the method of non zero gluon mass |^, however in this case we find that the expression 
for the power corrections involves two arbritary scales /(Ag^^/Q^) and J{Aq^jj/Q‘^) rather 
than a single gluon mass parameter A. The predictions for the z dependance will be different 
depending on which prescription we choose. The origin of these two different functions is 
that J arises from integrals which are collinear divergent and I from those that do not have 
this divergence. It is therefore not surprising to expect different functions appearing in this 
manner. This situation is not improved if we let the coupling in the various terms of Eq.(^) 
run in a different manner depending on whether they come from collinear or soft regions as 
z^l. 

There are two aspects of the above discussion that we would like to emphasize. (1) The 
predictions for the observable z dependence of the power corrections depend crucially on the 
argument of the running coupling. For the case of F^, if the coupling runs as as(k^) then the 
results disagree with the case 7 ^ 0 while the running as reproduces the results of 

the previous section ( see also [^). The argument of the running coupling can be clarihed 
in perturbation theory by two-loop calculations and has been studied for say the structure 
function F 2 in a number of papers as mentioned above, HI. (2) More importantly, we see 
that for the structure function F 2 the renormalon technique necessitates the introduction of 
two unknown scales(associated respectively with I and J) rather than just one in the case of 
A^ 7 ^ 0. Thus keeping in mind the transition to the non-perturbative case we should reserve 
for at least two independant rescaling functions. In this respect the picture is different from 
that in section 2. We will see a furthur indication of this in the next sections. 


4 Quark Mass as an Infrared Parameter. 

In an attempt to further quantify the infrared sensitivity in deep inelastic scattering phe¬ 
nomena, in this section, we consider the quark mass as an infrared parameter substituting 
say the gluon mass of section 2. It is convenient to discuss the role of the quark mass for the 


12 




moments of the structure functions, i.e. in the language of OPE. Moreover in this language, 
as discussed at the end of this section, it is most convenient to hnd the possible sources 
of infrared instability. In order to make the connection to earlier work m which also ex¬ 
ploited the OPE more transparent we choose to consider in this section, the moments of the 
antisymmetric structure function F^. 

The forward Compton amplitude has the standard operator product expansion: 

' 2 ^ ” 


- fll' 


= i (d^xexp{iqx)T{i^{x)il{Q)) = ^ 


-r 


■ - (35) 


{9tJ.fii9ufj,2Q 9tiiJ,iQuQ^i2 9vfi29.ii9.iii T 9fiv9fii9fi2)^2 ^^ivfiia9a9fi2^'^9fi3---9fin^fi^...fi„ 


where C* are the coefficient functions and the 9 are the operators. We will be interested in 
the 1/Q^ correction keeping the leading twist contribution as well. The calculation of section 


2 which reproduces the results of Ref. corresponds to a two-step procedure in evaluating 
the matrix elements of the operators 9 [|T^. Namely, one hrst reduces the operators of higher 
twist to those of the leading twist and then expresses the latter in terms of phenomelogical 
structure functions in a standard way. 

We will explain this procedure in detail on the example of the lowest moment of F^. The 
relevant operator is 0 , |^: 


'^fiv ^2^iiuaq9a f-^/3 + ^3 


(36) 


where 


Lp = 9199, 


9c = 9s9GlgFi9i59, 


- 2 




(37) 

(38) 

(39) 


Now, we assume that the matrix element of the operator 6^/3 containing the gluon held can be 
reduced to that of by using perturbation theory: i.e, we simply evaluate the contributions 
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from the diagrams of Fig.l with quark external states. (The open circles represent the 
insertion of the higher twist operator). In this way we get, 

Cp 4:0 s 


(P|r4,b> « ({P|i/3|P) - 


(40) 


where the factor A^) already mentioned in the Introduction is: 

/(e^,m^,A^) = ^ f dxM'^ln(M^), = — e^x(l — x) + + A^(l — a:). (41) 

Jo 

The evaluation of the l/Q'^ correction considered in section 2 |^, |T^, |T^ corresponds to 
/(0,0,A2) = -2(AVQ2)/nA2. 

Eq. ( 0 ) demonstrates an origin of the renormalon ambiguities in DIS. From the point of 
view of the general OPE (see Eq. (^6])) there is no reason whatsoever to treat the operators of 
the quark and gluon helds in an asymmetric way. The assumption that one can Erst integrate 
out the gluon line perturbatively and use a non-perturbative input on the matrix element 
of quark helds is arbitrary. This can be contrasted with the basic idea of renormalons[Q, Q 
when only one gluon line is made soft by means of an insertion of a large number of vacuum 
bubbles. Now we have both gluon and quark lines soft. It is only in this way that one can 
get 1/Q^ corrections. Thus the estimates of the 1/Q^ corrections in DIS can be understood 
only within a more general framework when one introduces infrared sensitive parameters in 
various ways. Then the quark mass is no worse a parameter than the gluon mass and one 
can look for m?lnrn? terms instead of terms. However for the case of the quark mass 

must be careful to take into account the purely kinematic higher twist contributions a la 
Nachtmanl^. To illustrate this point let us consider the second moment of the structure 
function F^. 

The relevant operator now is 0, |^: 


^^^i/aqqaq-yiLpry 


7 


■-S/ 37 ) 


(42) 


q- . 32g2 I6g2' 

where, the leading twist operator is also included and symmetrization and trace subtractions 
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are not explicitly shown. Furthur, 


= iqinD^q 
Bfh = gq[G%„, sB.,}+«“7„75g 


(43) 

(44) 

(45) 


Just as before, the matrix elements of the operators containing the gluon held can be reduced 
to that of by evaluating the contributions from Fig.l. The results for the two operators 
are: 


It \ Cl p ij,uafSQ.aQ.iy 

{Aq,) = —a.--- 


7?(p^m^A^) + /^“(p^m^A^) ^ 

(f ) 

, Cp^ f\p\m\X^) 

r) '^■S 2 

ZTT 


(Lq.,) (46) 

{h-,) (47) 


where, 



/f (p^, A^) = 9 / dx{l — x)M‘^lnM‘^ 

Jo 

(48) 


= —3 / dxx(—p‘^x(l — x) + m^x)lnM‘^ 

Jo 

(49) 


/^(p^, m^, A^) = —28 [ dxxM'^lnM'^ 

Jo 

(50) 

It is easy to verify that for 

= w? = 0 the results of ref. 13 are reproduced. 

For hnite 


quark masses, we see from Eq.(^6D and Eq,(^9D, that it appears in a different way than the 
gluon mass through the function f^. The origin of this function however is easily understood 
as a consequence of a kinematic effect due to hnite quark masses. To see this we hrst note 
that in general the coefficient function C{q) will depend on the running quark mass. 


m{Q) = m(Qo) + m{Qo)'ymln^ 

Vo 


(51) 


where, 7 ^ is the mass anomalous dimension which at the one loop order is: 


7m 


—o —as 
271 


(52) 
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The coefficient function itself may be expanded: 


C{q) = C°{q, m) + 771^(7° + .... 

m am 


(53) 


where the ellipses denote other higher order terms. The mass dependance of the coefficient 
functions can be inferred to any order in m/Q by considering the structure functions as a 
function of the Nachtman variable |]^ 




(54) 


Using this it is straightforward to reproduce the contribution /^. Thus we see that when the 
purely kinematic effects due to the quark masses are taken into account following a procedure 
analogous to that of Ref. [0, we get similar predictions for the infrared sensitivity keeping 
m as an infrared parameter or 7 ^ 0 , at least upto the second moment of T 3 . 

The formulation discussed in this section is the appropriate language for the inclusion of 
the anomalous dimensions of the operators governing the power corrections. A discussion of 
this, as we see next,brings into question the infrared safety of the relations like those in Eqs. 
( 0 ), ( 0 ) and (ID. 

Let us go back to the case of the hrst moment of T 3 . The leading twist contribution and 
the first power corrections are governed then by the operators L/j and respectively (see 
Eqs. ( p7t |38[)). The anomalous dimensions of the operator is obviously vanishing while 
that of the operator is []^: 

/ 32 , /" c: cN 

( 55 ) 

Since the anomalous dimensions of L^, and 6^ are different, Eq. (0) can at best be true 
only for a particular choice of . In other words, the relations like Eqs. (0) and in general, 
(0), ( 10) and are not infrared safe. If one accounts first for a splitting of the quark into 
a collinear quark and gluon and then for the emission of a soft gluon, the corresponding 
contribution to the power correction is not suppressed by q;s(Q^). This is the meaning 
of a nonvanishing anomalous dimension of Op. Note that in this respect the properties of 
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the power like corrections differ from the properties of the leading twist contribntions. In 
the latter case there are no new collinear divergences in higher loops if one considers the 
non-singlet DIS strnctnre fnnctions. This in tnrn means that the power corrections are 
governed by an independent strnctnre fnnction and its rednction to the leading twist one 
cannot be justihed in any known approximation ( for an earlier discnssion, see |P). As 
mentioned in section 3, an indication of this is already present in the renormalon chain 
method of obtaining the power corrections where these were fonnd to involve two different 
scales throngh the integrals / and J. 


5 Conclusions. 


The pictnre of predictions for the power-like corrections to deep inelastic scattering is remark¬ 
ably rich. From a pnrely phenomenological viewpoint, one deals not with a few nnmbers 
hxed by renormalons as is nsnally the case bnt rather with a few fnnctions which are x- 
dependences of the power corrections to varions strnctnre fnnctions []^ [^, |^. This makes 
comparison with the experimental data mnch more challenging. On the theoretical side, 
as is emphasized in this paper, there is an nnnsnal variety of possible choices of infrared 
sensitive parameters. We considered in detail the conseqnences of choosing the glnon and 
qnark masses as well as the position of the Landan pole in the rnnning conpling as infrared 
parameters. To snmmarize the resnlts, let ns mention some of the conclnsions made: 

(1) Keeping the glnon mass A 7 ^ 0 and defining the infrared sensitive contribntions as 
those non-analytic in A^, resnlts exactly in the same predictions as the dispersive approach 
to the conpling |n| and the renormalon chain in the large Nf limit m- The choice between 
these techniqnes is pnrely a matter of convenience and the case of A 7 ^ 0 looks most simple 
from the compntational point of view. 

(2) Parametrizing the power correction in terms of Aqcd brings ont two interesting 
featnres. The predictions based on A 7 ^ 0 are reprodnced if one assnmes that the conpling 
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runs as as{k\/{1 — z)). On the other hand, only by letting the coupling run as as{k\) do 
we get consistency with explicit two-loop calculations in case of the structure function F 2 . . 
This latter choice results in a different pattern of the power corrections to F 2 and to Fl than 
the case A 7 ^ 0. The power correcttions to Fl{x, are more important phenomenologically 
since Fl{x, Q"^) is proportional to in fhe leading-twist approximation. In fact here at 

the two loop level, one would encounter diagrams of the type shown in Fig. 2 which can be 
shown to have power corrections of the type at the partonic level. It is not clear that 

such contributions can be neglected even though they arise at a higher perturbative order. 
We also found that because of the collinear divergence , not one but two unknown scales 
are , in general, introduced in the renormalon chain method. One which is associated with 
h • cind the other with Ji ■ ■ This again is different from the case A 7 ^ 0 and is 

indicative of an inherent infrared instability. 

(3) The use of the quark mass to identify the infrared sensitive contribution results, 
generally speaking, in a different pattern of power-like corrections. The use of the hypothesis 
on non-perturbative enhancement of the infrared sensitive contributions Q] is crucial at this 
point. Namely one splits im?ln m terms into two pieces. One piece is a power-like corrections 
due to a finite mass of the target- the quark in our case. Such contributions may be treated 
in a manner similar to that in . The other piece is to be treated as a signal for an infrared 
sensitive contribution associated with low k^. Allowing for a non-perturbative enhancement 
of this piece brings the predictions in line with the case A 7 ^ 0, m = 0. 

(4) From the point of view of the OPE approach, the calculations of sections 2, and 
3 correspond to an asymmetrical treament of the effects due to the soft gluon and quark 
lines. One integrates over the gluon line perturbatively while the quark distribution (as 
manifested in the structure functions) are borrowed from experiment which implies the in¬ 
clusion of both perturbative and non-perturbative infrared effects. The dependence on the 
quark virtuality inherent to the reduction factor (^iD signals that the procedure cannot in 
fact be substantiated theoretically. However, the virtuality of the quark is not a convenient 
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parameter because it is manifestly not gauge invariant. To circumvent this we discussed the 
problem in the framework of the anomalous dimension of the operators governing the power 
corrections. Since these are different for the leading twist and the higher twist operators, we 
concluded that in general we are dealing with new infrared unsafe quantities at the level of 
the power corrections which cannot be cured by just introducing the same structure function 
as for the leading twist. Hence independant structure functions are required for the power 
corrections. Thus, experimental conhrmation of the one-loop results given in and 

reviewed in section 2 of the paper would in fact bring about a puzzle because of a successful 
perturbative evaluation of an infrared unsafe quantity. One would then be confronted with 
the challenge of formulating the approximation involved in more precise terms. 
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